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Abstract 

We do several simple calculations and measurements in an effort to 
gain understanding of global warming and the carbon cycle. Some 
conclusions are interesting: (i) There has been global warming since 
the end of the "little ice age" around 1700. There is no statistically 
significant evidence of acceleration of global warming since 1940. (ii) 
The increase of CO2 in the atmosphere, beginning around 1940, accu- 
rately tracks the burning of fossil fuels. Burning all of the remaining 
economically viable reserves of oil, gas and coal over the next 150 
years or so will approximately double the pre-industrial atmospheric 
concentration of CO2. The corresponding increase in the average tem- 
perature, due to the greenhouse effect, is quite uncertain: between 1.3 
and 4.8K. This increase of temperature is (partially?) offset by the 
increase of aerosols and deforestation, (iii) Ice core samples indicate 
that the pre-historic CO2 concentration and temperature are well cor- 
related. We conclude that changes in the temperatures of the oceans 
are probably the cause of the changes of pre-historic atmospheric CO2 
concentration, (iv) Data suggests that large volcanic explosions can 
trigger transitions from glacial to interglacial climates, (v) Most of 
the carbon fixed by photosynthesis in the Amazon basin returns to 
the atmosphere due to aerobic decay. 

1 Introduction 

We, two non-experts, present several "back-of-the-envelope" calculations and 
some simple measurements related to global warming and the carbon cycle. 
Our purpose is to understand which phenomena are important in determining 
the temperature of the Earth. 
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2 Data 



In this article we will use the following measured data. 0] The power of the 
radiation of the Sun per square perpendicular meter above the atmosphere 
is measured to be Ig U n — 1368W/m 2 . 1 This is the "solar constant". Since 
the surface of the Earth is 4 times the surface of a disk of the same radius, 
the average incident solar power per square meter of the Earth's surface 
is 1368/4 = 342W/m 2 . All fractions listed below refer to 342W/m 2 . The 
fraction of sun light that is reflected by the Earth is measured to be 0.31 
(0.21 by clouds, 0.06 by air including dust and water vapor, and 0.04 by the 
ground). This is the "albedo" of the Earth. The remaining fraction (0.69) is 
absorbed by the Earth (168W/m 2 by the surface, 48W/m 2 by the troposphere 
(water vapor and aerosols), 10W/m 2 by the stratosphere, and 10W/m 2 by 
clouds). So, the net solar power absorbed by the Earth is 342 • 0.69 = 
236W/m 2 for average cloud coverage. 

The infrared radiation of the Earth to space measured by satellites is 
235 W/m 2 for average cloud coverage (88W/m 2 by clouds, 126W/m 2 by water 
vapor and CO2, and 21 W/m 2 by the surface through the infrared atmospheric 
window). So, the incoming and outgoing energies balance to high accuracy 
and determine the temperature of the biosphere. 

The average temperature of the surface of the Earth is 288K= 15°C. The 
average temperature of the atmosphere is about 250K= — 23°C. 

In cloud free conditions the average solar power absorbed is 286W/m 2 , 
while the average emitted terrestrial radiation is 266W/m 2 . The correspond- 
ing numbers for average cloud cover are 236W/m 2 and 235 W/m 2 as men- 
tioned before. Note that decreasing the cloud cover has a net heating effect. 

The Earth radiation spectrum measured from a satellite (Figure P) shows 
that the infrared windows (A = 8 to 9/zm and 10 to 13/xm) radiate as a 
blackbody at the surface temperature, the CO2 absorption band (A = 14 to 
16/im) radiates at 215K (corresponding to the altitude of the tropopause at 
12km), and the H 2 bands (A < 8/im and A > 17/xm) radiate at ~ 260K 
(corresponding to an altitude of 5km). 

lr The power per unit area radiated by a black body at absolute temperature T is erT 4 , 
with a = 5.67 • l(T 8 Wm- 2 K- 4 . The power radiated by the Sun is P = ^R 2 Sun aT^ un . 
The power per unit perpendicular area received by Earth (outside of the atmosphere) is 
1368W/m 2 = -P/(47ri?g un Earth ). From these equations we obtain the effective temperature 
of sun light: T Sun = 5780K. 
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Figure 1: Radiation spectrum of the Earth as observed by a satellite over 
Africa. [3] Wavenumber is 1/A. Blackbody spectra at several temperatures 
are also shown. Multiply the numbers on the vertical axis by 7r to obtain 
the radiation in 10 _3 W m" 2 (cm -1 ) -1 (with m 2 now referring to the surface 
of the Earth instead of the satellite antenna). 
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3 A simple model 



Let us begin with a simple model. From observations we know that, to high 
accuracy, there is an equilibrium between incident solar power and radiated 
infrared power. Therefore 

-^vlsun = tlRCT^, (1) 

where Is w = 1368W/m 2 is the solar constant, e v is the fraction of incident 
sun light power that is absorbed by the Earth, i.e. it is the "emissivity in 
the visible", em is the emissivity in the infrared, a = 5.67 x 10~ 8 Wm -2 K -4 
is the Stefan-Boltzmann constant, and T is the temperature of the surface 
of the Earth in degrees Kelvin. The factor | was explained in Section 2. 

We note that if e v = em the equilibrium temperature is 279K, quite close 
to the observed mean temperature of 288K. 

From the data presented in Section 2 we conclude that e v ~ 0.69 for 
the Earth as a whole, e v ~ 0.55 for areas covered by clouds (it varies from 
0.6 for cirrus to 0.1 for cumulonimbus), e v ~ 0.9 for cloudless ground, and 
e v ~ 0.2 for snow. Because snow is mostly at high latitude we sometimes 
replace e v ~ 0.1 for snow. 

From the data presented in Section 2 we estimate that em « 0.605 for 
the Earth as a whole, em ~ 0.56 for areas covered by clouds, em ~ 0.69 for 
cloudless ground, and em ~ 0.3 for snow. 

These "effective" emissivities are valid for the model of this Section, i.e. 
an Earth surface characterized by a single temperature, and includes the 
atmosphere with its greenhouse gases. The model is too crude to account for 
all observations, so these effective emissivities are approximate. A summary 
of effective emissivities is presented in Table [0 In the last column of the 
table we show the equilibrium temperature for 100% coverage of black body, 
cloudless ground, clouds or snow. The last row shows the world average. 

4 Two box model 

Figure El shows two boxes: one for the upper atmosphere at temperature 
T2, and one for the surface and lower atmosphere at temperature T.jl] The 
surface radiates ~ <rT 4 . A fraction 1 — e of this radiation gets through the 
upper atmosphere (this is the infrared window), and a fraction e is absorbed 
by greenhouse gases and aerosols in the atmosphere. The upper atmosphere 
radiates ecrT 4 back to the lower atmosphere, and eoT^ out to space. The 
total radiation out to space is 

aT A - ea (T 4 - T 4 ) . (2) 
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1.00 


279K 


cloudless ground 


« 0.69 


0.90 


298K 


clouds 


« 0.56 


0.55 


278K 


snow 


« 0.30 


0.10 (0.20) 


212K (252K) 


whole Earth 


0.605 


0.69 


288K 



Table 1: Effective emissivities for infrared and visible radiation. The last 
column shows the temperature corresponding to a 100% cover of black body, 
cloudless ground, cloud or snow. The last row corresponds to 50% cloud 
cover, 5% snow cover, leaving 45% clear ground. The measured average 
temperature is 288K. 



The last term reduces the infrared emission of the Earth (at a given T) 
because the upper atmosphere is colder than the surface. This last term is 
the greenhouse effect. The effective emissivity is 



T 



(3) 



We estimate e ~ 0.813 (see next section). For T = 288K we obtain T 2 ~ 
244K. Note that heating of the upper atmosphere reduces the greenhouse 
effect. Adding greenhouse gases to the atmosphere has two opposing effects 
on T: e increases and T 2 increases. 



5 Three box model 

The three box model is shown in Figure El The boxes are the surface of 
the Earth at temperature T, the lower atmosphere at temperature T\, and 
the upper atmosphere at temperature T 2 . The total radiation to space is 
equal to the total solar radiation absorbed (235W/m 2 ), which we assume 
independent of the concentration of CO2. The solar power absorbed per unit 
area is 168W/m 2 by the surface, and 67W/m 2 by the atmosphere. Thermals, 
evaporation and transpiration are estimated to transport 102W/m 2 .[lJ The 
atmosphere absorbs radiation at frequencies near vibration resonances of the 
greenhouse gases, mainly H 2 0, CO2 and ozone. The fraction of infrared 
radiation absorbed is e. The atmosphere also radiates at these frequencies. 
The emissivity of the surface in the infrared is taken to be 1. The sum of 
powers entering each box is zero. The greenhouse effect is due to the low 
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Figure 2: Two-box model of the greenhouse effect. ^ 



temperature of the upper atmosphere, so that emission to space is reduced 
at frequencies near resonances as shown in Figure H 

In Table |5] we show the temperatures as a function of the emissivity of 
the atmosphere in the infrared (e). The absorbed sun power and power of 
thermals, evaporation and transpiration are assumed to be constant (inde- 
pendent of the concentration of CO2), so the model breaks down for e < 0.5 
(with no atmosphere we should obtain T = 253. 72K for 235 W/m 2 black-body 
radiation) . 

It has been estimated that doubling the concentration of CO2 from (290ppm 
to 580ppm) decreases the radiation to space by 3.75 ± 0.25W/m 2 1J (before 
temperatures are allowed to adjust) which corresponds to increasing e from 
0.813 to 0.833. Note, in Table El that doubling the concentration of CO2 
raises the surface temperature T by 2. IK. Do we trust this result? No! Most 
of the radiation to space is from the upper atmosphere. The temperature 
of the upper atmosphere may change with CO2 concentration due to energy 
fluxes or absorptions not considered in this simple model. Also the solar 
inputs and power of thermals, evaporation and transpiration surely depend 
directly or indirectly on the CO2 concentration (by changing T and 7\). 
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Figure 3: Three-box model of the greenhouse effect. The numbers are in 
units of W/m 2 . Dotted lines correspond to the infrared window, while dashed 
lines correspond to frequencies absorbed or emitted by greenhouse gases in 
the atmosphere (mainly H 2 0, CO2 and ozone). Thermals, evaporation and 
transpiration are estimated to transport 102 W/m 2 . 
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Table 2: Temperatures of the surface (T), lower atmosphere (7i), and upper 
atmosphere (T 2 ) in the three box model as a function of the emissivity of the 
atmosphere in the infrared (e). The absorbed Sun power and power of ther- 
mals, evaporation and transpiration are assumed to be constant (independent 
of the concentration of C0 2 ). 

6 Global warming from Earth emission spec- 
tra 

Let us try a different approach that relies only on the Earth emission spectra 
shown in Figure ^ We estimate the increase in surface temperature due to 
a doubling of the concentration of C0 2 (from 290ppm to 580ppm). We take 
the average Earth temperature to be 288K (instead of the 320K of Africa 
shown in Figure Q). 

C0 2 has two absorption bands in the tail of the solar spectrum, which 
absorb ~ 0.8W/m 2 . We will neglect this absorption (a correction could be 
applied later). Increasing the concentration of C0 2 will widen the absorption 
resonance seen in Figure ^ How will the spectra respond? If the Earth 
albedo remains constant, then the area below the spectra in Figure ^ will 
remain constant. Then the spectra will respond by increasing (or in general, 
modifying) the surface temperature, and/or the water emission temperature 
(Th 2 o ~ 260K), and/or the C0 2 emission temperature (Tco 2 ~ 215K). Table 
El suggests that T\ and T 2 vary less than the surface temperature T. So, in 
this section, without much justification, we will assume that (i) Tjj 2 o and 
Tqo 2 remain constant, and (ii) the Earth albedo remains constant. So, in 
this approximation, the only response to the widening of the C0 2 absorption 
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resonance is to increase the surface temperature to conserve the area under 
the spectrum in Figure ^ 

Let us estimate the widening of the CO2 resonance using only the data 
in Figure E This widening in W/m 2 (before any temperature has a time to 
react) is called "radiative forcing". We label points of transmittance equal 
to 1, 0.7, 0.5, 0.25, 0.06, 0.004 and 0.0 along the side slopes of the absorption 
resonance, and then shift them downward (at constant wave number). 2 The 
result is a radiative forcing « 4W/m 2 . The corresponding change in surface 
temperature, in responce to a doubling of the concentration of CO2, with the 
assumptions discussed above, is ~ 2.5K. 

Another estimate would let T and Th 2 o vary by the same amount, while 
Tqo 2 remains constant. The corresponding warming is ~ 1.3K. 

7 Detailed models 

The results presented in this section were obtained from £Q. The surface and 
troposphere are tightly coupled (by non-radiative heat exchanges), while the 
coupling between the troposphere and stratosphere is relatively weak. For 
this reason changes in the surface temperature are largely determined by 
changes in the net (incoming minus outgoing) radiation at the tropopause 
(at an altitude between 10 and 20 km). This change in net radiation at 
the tropopause per unit surface area (before any temperatures are allowed to 
change), AR, is called "radiative forcing" . The change in surface temperature 
is expressed as AT = AR/X, where A is the "radiative damping". If the 
Earth were a blackbody (radiating the same amount as the Earth), then 
A = 3.7Wm~ 2 K _1 . However due to various feedback effects, it is estimated, 
using detailed models, that A = 2.0 ± 0.5Wm" 2 K _1 .jT] So the net feedbacks 
are positive. 3 

The change in radiative forcing and surface temperature, due to changes 
in composition of the atmosphere since the pre-industrial period until 1995, 
is estimated by detailed modelspQ to be AR = 2.6 ± 0.6Wm~ 2 and AT = 
1.3 ± 0.44K due to all greenhouse gases (about half of this is due to Pco 2 
increasing from 290ppm to 350ppm), and AR = — 1.9 ± l.lWm" 2 and 
AT = —1.0 ± 0.6K due to all aerosols (we have added errors in quadrature). 
Deforestation changes the albedo of the Earth, causing a change of tempera- 
ture ~ — 0.4K.[1] So we do not know wether the increase of greenhouse gases, 
aerosols and deforestation has a net heating or cooling effect. 

2 To find where to place the points and by how much they shift downward, we simulated 
a ten-layer atmosphere model on a spreadsheet. 

3 Why A drops can be understood from our discussion in Section 6. 



9 



Let us now consider a doubling of the pre-industrial concentration of CO2 
(from 290ppm to 580ppm). The change in radiative forcing is AR = 3.75 ± 
0.25Wm _2 pQ and the corresponding change of surface temperature is AT = 
1.9 ± 0.5K. This result, obtained from detailed models^] is in agreement 
with paleo-climate and historical data (and with our back-of-the-envelope 
estimates of Sections 5 and 6). In comparison, the change in temperature 
since the "little-ice-age" (1700) to the present (1995) has been pa 1.2K.j3] 
It is worth mentioning that various detailed models obtain AT in the range 
from 1.3K to 4.8K for doubling of The range of results is large so this 

is a difficult problem that does not seem to converge: as the model becomes 
more complex, the uncertainties appear to increase! 

8 Ice ages 

For given emissivities, the Earth temperature is stable because an increase in 
temperature causes an increase in the infrared radiation, which in turn brings 
the temperature back to its equilibrium value. However, the equilibrium 
temperature depends on the emissivities. For example, a world-wide snow 
storm could result in a shift of the equilibrium temperature from 288K to as 
low as 252K, see Table ^ which would prevent snow from melting. In this 
way an ice age might begin. 

The mean annual temperature averaged over the Earth fluctuates from 
year to year with a standard deviation of about pa 0.22°C (during the 20'th 
century). Inter-glacial periods last typically pa 50 thousand years. With 
50% probability we expect a fluctuation of at least 4.2 standard deviations 
in 50000 tries. This corresponds to pa 1°C. Is such a world-wide temperature 
fluctuation sufficient to throw us into an ice age? 

To get out of an ice age, an event is needed to change the emissivities, 
perhaps volcanic eruptions (or meteorite impacts) that cover snow with ash. 
Table El is a list of the largest known volcanic eruptions in the last million 
years. Each one of them (except perhaps Taupo in isolated and tropical 
New Zealand) corresponds (within errors) to a transition from a glacial to an 
interglacial period, see Figure EJ It is interesting to note that an explosion 
that ejects 2 x 10 15 kg corresponds to an average of 4kg-m -2 over the entire 
Earth, so that such an event could have a sizable effect on the albedo. If 
the volcanic origin of the transitions is correct, we should find evidence for 
other large eruptions at approximately 770, 430, 340, 270 and 140 thousand 
years ago. Note however that volcanic explosions have a short-term cooling 
effect. |U 

On 4 November 2002 the volcano El Reventador in Ecuador exploded and 
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Toba (Indonesia) 


6.9 


74 


Yellowstone (Wyoming, USA) 


2.2 


600 


Porsea (Toba, Indonesia) 


2.0 


790 


Taupo (New Zealand) 


1.3 


26.5 


Long Valley (California, USA) 


1.2 


700 



Table 3: Large explosive volcanic eruptions in the last million years. [H| 

covered Cayambe with ash. Two years later we see a barren Cayambe with 
large patches of glacier gone. Similarly, volcanic ash from Tungurahua has 
been falling on Chimborazo for the last few years. Chimborazo also looks 
barren compared with what it used to be. 



9 Width of the 15/rni line of C0 2 . 

The CO2 molecule is linear, and has three normal modes of oscillation. They 
have wavelengths of 15/zm, 7.46//m and 4.26//m. The mode with wavelength 
7.46/mi has no dipole moment and does not couple to electromagnetic radia- 
tion. The mode with wavelength 4.26/im barely overlaps the black-body radi- 
ation spectrum of the Earth. So the important spectral line has a wavelength 
A = 15/im. We calculate the half-width of this line (neglecting molecule col- 
lisions) at the points of half power to be AA/A ~ 10~ 13 . 

Let us consider molecule collisions. The calculated mean free path at a 
pressure of one atmosphere is I ~ 1/zm, the rms thermal velocity of CO2 is 
v = 400m/s at 300K, and the mean time between collisions is r pa 2.5ns. 
If the time between collisions is t, then the half-width of the spectral line 
between points of half power is AA/A ~ 0.45A/(ct). For t = 2.5ns the result is 
AA/A pa 10~ 5 . The probability that t is in the interval dt is r 1 exp [— t/r] ■ dt . 
So the spectral line has very long tails that absorb electromagnetic radiation. 
These long tails explain the wide absorption band of CO2 seen by satellites. 

We have calculated the half-width AA of the window corresponding to 
the A = 15/^m line of CO2 at which the power per unit area getting through 
the atmosphere drops by a factor e~ 2 : 



2Pco 2 Lrio 
kTmr 



(4) 



where rj = 1207T Ohm is the impedance of free space, Pco 2 is the partial 
pressure of CO2 at sea level, L = kT / (Mg) = 5550m is the altitude at which 
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Figure 4: Ice-core temperature data of the last 900 thousand years com- 
pared with the largest known volcanic explosions in that period (arrows) [S] . 
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Table 4: The columns contain the ocean pH; the number of kilograms of 
carbon per cubic meter of ocean at 0°C and 25°C in the form C0 2 and 
-ffCOjjZj; and the corresponding total carbon in the upper 1000m of the 
oceans. The calculations correspond to 300ppm of C0 2 in the atmosphere. In 
comparison, the atmosphere with 300ppm of C0 2 contains 430Gt of carbon. 
Oceans have a pH between 7.4 and 8.3. [7 

the partial pressure of C0 2 drops by a factor e _1 , c is the velocity of light, m 
is the reduced mass of a carbon atom oscillating against two oxygen atoms, 
T is the absolute temperature, r is the mean time between collisions of C0 2 
molecules, and q is the effective charge of the carbon in C0 2 . Substituting 
numerical values we obtain AA ~ 0.6/im for 300ppm C0 2 , in reasonable 
agreement with 1.0/im measured by satellites. So we understand the long 
tails of the spectral line of C0 2 . Note that doubling Pco 2 increases AA by a 
factor V2- 4 

10 Equilibrium of CO2 between atmosphere 
and oceans 

The chemistry of carbon in the oceans is beautifully described in We are 
interested in these two equations of chemical equilibrium: 

[C0 2 ] = K H P C02 , (5) 

which describes the solution of C0 2 in sea water, and 

[H + ][HCO^] = K al [C0 2 ] (6) 

that describes the ionization of hydrated carbon dioxide. Pco 2 is the partial 
pressure of C0 2 in the atmosphere (in atmospheres), and [...] are concen- 
trations in sea water (in mol per liter). The equilibrium constants have 

increasing AA by a factor results in a radiative forcing greater than our estimate 
of Section 6, so the calculation should be taken with a grain of salt when applied so far 
out on the tails of the resonance. 
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these values at 0°C (25°C): -\ogK H = 1.16 (1.54), and -\ogK al = 6.07 
(5.86) for 0.7 ionic strength. [7] The calculations shown in Table 0] assume 
P C02 = 0.0003atm (z.e.300ppm). 

Assuming that the total carbon in the air and the top 1000m of the oceans 
is constant, and that 300ppm corresponds to 10°C and pH = 7.4 (8.3), we 
obtain that a 10°C increase in the sea water temperature produces an increase 
of CO2 in the atmosphere of 40ppm (48ppm). In comparison, from Antarctic 
and Greenland ice cores, it is found that 9.5°C changes in temperature have 
been accompanied by ~ 75ppm change in Pco 2 - 

In conclusion, we qualitatively understand the coupling of the pre-historic 
temperature with the atmospheric concentration of CO<i- an increase of ocean 
temperature reduces its solubility to CO2, causing an increase of the atmo- 
spheric concentration of COi- The converse is not true: an increase of 75ppm 
CO2 would increase the temperature, due to the greenhouse effect, by only 
« 0.8°C. 

Let us now estimate the effect of acid rain due to a large volcanic ex- 
plosion. Assume that 2 x 10 15 kg are ejected. Assume that 0.1% of this is 
sulfur (in any chemical form) that ends up as sulphuric acid in the oceans. 
The concentration of H2SO4 in the top 1000m of the oceans would be 
1.8 • 10~ 4 mmol/liter. We measured the pH of ocean water off the coast of 
Ecuador as a function of H2SO4 concentration: 7.3, 7.0, 6.3, and 5.5 for 0, 
0.0084, 0.0211 and 0.0422mmol/liter respectively. 5 So the change of pH is 
negligible (~ 0.006), and the corresponding release of C0 2 by the oceans is 
also negligible. 

11 A comment on the carbon cycle 

It is observed that the concentration of CO2 has increased linearly from 
325ppm in 1970 to 375ppm in 2004. [HJ If this increase continues at constant 
rate it will take about ~ 140 years to reach twice the pre-industrial concen- 
tration of 290ppm. 

It is observed that the buildup of atmospheric Pco 2 accurately follows 
the human emissions of CO2 (due to burning of fossil fuels and the clearing 
of land for agricultural and urban use).pQ The factor of proportionality is 
observed to be 65%. The remaining 35% of emitted carbon apparently ends 
up in the surface (mixed) layer of the oceans and in bio-mass. No saturation 
is observed so far, so the time constant with which carbon reaches the deep 

5 Measurement done several days later on a sample of water taken from the beach at 
Atacames. A measurement done locally in Mompiche yielded pH = 7.7. 
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oceans is in excess of 50 years (it is estimated that there are several time 
constants, some as long as 2000 years). [T] 

There are enough world reserves of fossil fuels (gas, oil and coal) to last 
~ 150 years at the present rate of consumption (6.5Gt C per year). Burning 
these reserves would increase CO2 to three times the pre- industrial concentra- 
tion assuming 65% of the emitted C0 2 ends up in the atmosphere. However, 
on these longer time scales it is estimated that 30% (instead of 65%) of the 
emitted CO2 will remain in the atmosphere. pQ Thus we estimate that burn- 
ing all of the remaining reserves of fossil fuels in the next 150 years or so, 
will result in twice the pre-industrial CO2 concentration. 

12 Time constant to reach CO2 equilibrium 
with the oceans. 

The solubility of CO2 in the oceans decreases with increasing temperature. 
Consider a drop in the temperature of the oceans. With what time constant 
does the atmospheric concentration of CO2 drop due to the increased ab- 
sorption by the oceans? Here we assume that every CO2 molecule colliding 
with the ocean surface "sticks" to the surface. The number of CO2 molecules 
colliding with a square meter of ocean surface per second is[5| 

R = = 7 ■ lO^m-V 1 . (7) 

V 2vrM kT v ; 

The number of CO2 molecules in the atmosphere per square meter is Pao 2 /{Mg) = 
4- 10 25 m~ 2 . So the time constant for absorption of CO2 by very cold oceans is 
~ 90s (we have taken account that 0.7 of the surface of the Earth is covered 
by oceans). So equilibrium of CO2 and the oceans is very short if we neglect 
diffusion of CO2 through the atmosphere. 

Let us now consider diffusion. The partial pressure of CO2 decreases 
with altitude exponentially with a characteristic height L = kT j (Mg) = 
5550m. How long does it take a molecule of CO2 to diffuse 5550m through 
the atmosphere? The result is t = L 2 /(vl) ~ 2400 years. This result assumes 
no convection. Convection will reduce this time considerably. In fact, the 
flux of carbon from the atmosphere to the oceans and vice versa is enough 
to replace all the atmospheric carbon in 6 years. [1J So, on time scales of 
interest, the partial pressure of CO2 in the atmosphere is proportional to 
the concentration of CO2 in the surface water of the oceans: they are in 
equilibrium. 
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13 The carbon cycle in the Amazon rain for- 
est 



In the process of photosynthesis, a plant takes C0 2 from the air, breaks it 
up using solar energy, and releases 2 back to the atmosphere. The carbon 
attaches to water to form carbohydrates (CH 2 0) n , which is what plants are 
made of. So, for every molecule of oxygen released to the atmosphere, there 
remains one atom of carbon in the plant. What happens to the carbon when 
the plant dies? There are three alternatives: the carbon can (i) return to the 
atmosphere, (ii) accumulate on the ground, or (iii) end up in the oceans. 

The photosynthesis process can occur in both directions. Examples of 
photosynthesis in reverse direction are combustion when we burn wood, res- 
piration, and decomposition by aerobic bacteria. In these processes the car- 
bon of the plant becomes attached to oxygen from the air, and is returned to 
the atmosphere as CO2. Note that no net oxygen is produced: the oxygen 
released to the atmosphere during the plant growth is consumed during the 
plant aerobic decay. 

The carbon can end up in the soil, either as partially decomposed or- 
ganic matter (humus), as hydrocarbons or as coal. If the plant dies in a 
medium that lacks oxygen, such as in oceans, lakes, rivers, swamps or under 
volcanic ash, the anaerobic bacteria have a chance to decompose the plant 
(without the competition from aerobic bacteria). The anaerobic bacteria 
produce methane ("swamp gas"), carbonic acids, bicarbonates, and humic 
acids. The methane can escape to the atmosphere, or under enough pressure 
and temperature, can polymerize into hydrocarbons or can be crushed into 
coal. 

The carbon can be washed down the rivers, in the form of carbonic acids, 
bicarbonates, and humic acids which are soluble in water. The rivers in the 
Amazon basin are dark due to the humic acids. 

The Amazon rain forest has humus only in the first few centimeters of 
soil. So we neglect alternative (ii). To quantify alternative (iii), we took sam- 
ples of water of several tributaries of the Amazon river (Napo Alto, Napo 
Bajo, San Francisco and Tiputini) and measured the concentration of bi- 
carbonate (0.71 mol/m 3 ) and carbonic acid (0.11 mol/m 3 ). Assuming that 
these concentrations are typical of the Amazon river, and multiplying by the 
yearly discharge of the river, we obtain the carbon discharged to the Atlantic 
ocean by the Amazon river: 50 million (metric) tonnes per year. In total, 
the Amazon basin fixes about 10 billion brute tonnes of carbon per year. 
So, about 99.5% of this carbon returns to the atmosphere due to the aerobic 
decomposition of organic matter, while only about 0.5% is washed down the 
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rivers into the Atlantic ocean. 

It is interesting to compare these numbers with the world energy-related 
release of carbon dioxide from the consumption of oil, gas and coal: 24 billion 
tonnes in 2004. This corresponds to 6.5 billion tonnes of carbon per year. 
So the Amazon basin can transfer, from the atmosphere to the oceans, one 
year's worth of carbon dioxide emission in 130 years! 

14 Conclusions 

With some degree of confidence, we arrive at the following conclusions. 

1. Ever since we have accurate measurements of the temperature of the 
atmosphere, i.e. since 1702, we observe a global increase of tempera- 
ture of about 0.4°C per century. This warming does not accelerate in 
the second half of the XX century. In fact, we see no statistically signif- 
icant global warming since 1940. [2J|3] We observe no correlation of the 
global temperature with the consumption of oil, coal and gas, or with 
population. The atmospheric concentrations of CO2 begins to increase 
around 1940 [2]| when oil consumption takes off, yet we observe no corre- 
sponding increase of the slope of global warming. El ED] Some studies 
suggest that this is due to a balance between the heating effect of C0 2 
and the cooling effects of aerosols and deforestation. pill] From the data 
it is concluded that the accumulated effect of humankind on the global 
temperature until 1990 is not statistically significant: 0.0±0.1°C.|H] So, 
the observed global warming since 1702 is part of the natural variability 
of the climate as we pull out of the "little ice age" three centuries ago. 

2. The pre-industrial concentration of CO2 in the atmosphere was ~ 290 
ppm in 1800. It has since increased to 375 ppm in 2004. 8 There is 
convincing evidence that this increase is due to burning of fossil fuels 
(and forrests, and the manufacture of cement). 

3. Detailed models estimate that the changes of atmospheric composition 
since pre-industrial times up to 1995 produce a change of temperature 
of 1.3 ± 0.44K due to all greenhouse gases (about half of this is from 
CO2), — 1.0±0.6K due to aerosolsjl], and ~ — 0.4K due to the change of 
albedo caused by deforestation So we do not know if human activity 
heats or cools the Earth. As indicated above, no significant net change 
has been observed. 

4. The increase of the concentration of CO2 in the atmosphere accurately 
tracks the burning of fossil fuels (and forrests). Two thirds of the 
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carbon burned ends up in the atmosphere as CO2, and the remaining 
third ends up dissolved in the surface layer of the oceans (the time 
constant for this dissolution is about 6 years). The fraction of excess 
carbon in the atmosphere should decrease to one third in about 50 
years as more carbon is incorporated into the oceans (as CO2, HCO3 , 
CO^~, CaCOs), and in the land and marine biota. The last third of 
the excess C0 2 in the atmosphere is expected to return to the deep 
oceans and become buried as elemental carbon and as carbonates with 
several time constants ranging from hundreds to thousands of years. P] 

5. Burning all of the remaining economically viable reserves of oil, gas 
and coal over the next 150 years or so will approximately double the 
pre- industrial atmospheric concentration of C02- The global warming 
due to a doubling of the concentration of CO2 is expected to increase 
the average surface air temperature by 1.3K to 4.8K.pQ The increase of 
temperature is expected to be higher than average at higher latitudes 
(and lower than average at low latitudes). The heating effect of CO2 is 
obvious when we look at the Earth emission spectra shown in Figure H 
The range of results is large, so global warming is a difficult problem 
that does not seem to converge: as the models become more complex, 
the uncertainties appear to increase! 

6. Ice core samples indicate that CO2 concentration and temperature are 
well correlated: a change of temperature by 9.5°C corresponds to a 
change in CO2 concentration from 195ppm to 270ppm. 2] Does the 
change of CO2 concentration cause the change of temperature (by the 
greenhouse effect), or does the change of temperature cause the change 
of CO2 concentration (due to the temperature dependence of the sol- 
ubility of carbon in the oceans)? Due to the greenhouse effect, an 
increase of CO2 concentration from 195ppm to 270ppm would cause an 
increase in temperature of ~ 0.8°C, much too small to account for the 
observed change in temperature. Conversely, an increase of 9.5°C of the 
temperature of the upper 1000m of the oceans would increase the at- 
mospheric concentration of CO2 by m 44ppm (depending on the ocean 
pH). So it is plausible that the changes of temperature of the oceans 
caused the pre-historic changes of atmospheric CO2 concentration. 

7. Data suggests that large volcanic explosions can trigger transitions 
from glacial to interglacial climates, by covering ice with volcanic ash, 
thereby changing the equilibrium temperature of the Earth as seen in 
Table □ 
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8. Photosynthesis in the Amazon basin fixes ~ 10 billion tonnes of carbon 
per year. However, aerobic decay releases almost all of this carbon back 
to the atmosphere as CO<i- The balance, ~ 50 million tonnes of carbon 
per year, is discharged to the Atlantic Ocean in the form of carbonic 
acids, bicarbonates, and humic acids. 
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